On universality of algebraic decays in Hamiltonian systems 
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Hamiltonian systems with a mixed phase space typically exhibit an algebraic decay of correlations 
and of Poincare recurrences, with numerical experiments over finite times showing system-dependent 
power-law exponents. We conjecture the existence of a universal asymptotic decay based on results 
for a Markov tree model with random scaling factors for the transition probabilities. Numerical 
simulations for different Hamiltonian systems support this conjecture and permit the determination 
of the universal exponent. 
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The phase space of Hamiltonian systems with two 
degrees of freedom generally shows an intricate mix- 
ture of chaotic and regular structures. Regular regions 
consisting of 'islands' of quasi-periodic motion (KAM 
tori) appear hierarchically interspersed in a chaotic sea. 
The character of the motion in the irregular component 
is dominated by stickiness of trajectories close to the 
boundary circles separating regular from irregular re- 
gions. Stickiness crucially affects statistical quantities, 
like correlations, the distribution of Poincare recurrences 
or anomalous diffusion, which show related algebraic 
behaviors. Finite-time numerical experiments revealed 
system-dependent values for the power-law exponent z 
of the decay of the distribution of Poincare recurrences, 
P{t) ^ t^^, that is the probability to return to a given re- 
gion after a time larger than t. Great effort was devoted 
over the last two decades to understand if there is an 
asymptotic universal decay for P{t) and how its temporal 
behavior is related to the hierarchicalphase-space struc- 
ture near regular boundaries Q, @, i, i, d, @, 0, S, i, [TO] ■ 
In particular, renormalization techniques were used to 
derive scaling relations for the self-similar structure close 
to the critical circle with golden mean frequency [^, 0| 
or for trapping in exact self-similar island-around-island 
structures Q. More recently it was demonstrated 
that both scenarios are usually present in an Hamiltonian 
system, thus strengthening the validity of a Markov tree 
approach 0. While the self-similar Markov tree model 
succeeds to reproduce a power-law decay, it is based on 
the approximation of exact self-similarity 0, [Hj • It re- 
mains an open question how deviations from exact scaling 
affect the algebraic decay and the issue of universality. 




In this paper we study the role of random scaling 
factors in a Markov tree model. Initially, different 
realizations of the randomness give rise to quite different 
power-law behaviors (see Fig. [1]). In contrast, we show 
that asymptotically the exponent of the algebraic decay 
is the same for all realizations. We argue that this is rep- 
resentative for Hamiltonian systems and present strong 



FIG. 1: (Color online) (a) Survival probability Pi{t) for 100 
realizations of random scaling factors (thin lines) and their 
logarithmic average (thick line). Inset: binary tree with root 
s and absorbing site Ds. (b) Rescaled survival probabilities 
Pi(t)t^' with zi from Eq. (|lUp . showing the same asymptotic 
power-law exponents, (c) Variance cr^(t) of ln[Pi(t) t^'] show- 
ing saturation. 



numerical support. As a consequence the average of 
P(t) over different Hamiltonian systems is a meaningful 
quantity and allows for an efficient determination of the 
asymptotic universal decay. We numerically extract a 
universal exponent z ~ 1.57 for Poincare recurrences in 
Hamiltonian systems. 

We consider a Markov model on a binary tree. For 
an arbitrary site s, we call Ds the site obtained by go- 
ing up one level in the tree towards the root site s (see 
Fig. [T^, inset). The dynamics on the tree is defined by 
the probabilities to move in one time step to one of the 
neighboring sites. As we are interested in the survival 
probability we put Ds as an absorbing site and we fix a 
non-zero probability of leaving the tree pg^Ds UM ■ The 
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transition probabilities for all sites s below s are then 
chosen according to 



PDs^s/Ps->Ds = As 
Ps^Ds/PDs^DDs = Bs 



(1) 



where the scaling factors Ag and Bs are positive ran- 
dom variables with the constraints (i) (A) < 1/2 and (ii) 
B < 1. The constraint (ii) ensures that the transition 
probabilities decrease moving away from the root s. The 
constraint (i) allows to define an initial probability den- 
sity ps{t = 0), with normalization ^sP^i^) = 1; that 
would be invariant for the closed system {pg^Ds = 0) by 
satisfying detailed balance ps{0)ps^Ds = PDsiO)PDs^s- 
Due to the non-zero transition probability ps^Ds of leav- 
ing the tree towards the absorbing site Ds there is a de- 
caying survival probability 



Piit) :=E^-W' 



(2) 



of being still inside the tree at time t, with the subscript I 
standing for the "invariant" initial condition. 

In order to quantify the decay of the survival proba- 
bility we make use of the following observation: In the 
sum, Eq. ([2]), the density Ps{t) at a site s gives its most 
important contribution while it is decaying exponentially, 
with a rate that is dominated by the probability p^^Ds to 
escape in upward direction. At later times, even though 
its decay becomes algebraically slow, its contribution can 
be neglected [l^l- The survival probability can thus be 
approximated by a sum of decaying densities [3] 



Pl(t)«^p.(0)£ 



(3) 



Numerically, we find this to be a good approximation, in 
the sense that the ratio with the exact decay of Eq. ([2|) 
is found to be asymptotically constant. 

Power-law asymptotics can be conveniently extracted 
using singularity analysis of the Mellin transform. An 
algebraic decay P{t) ^ t~" is transformed into a simple 
pole of its Mellin transform P*{z) = J dtf"-^ P{t) at 
z ~ a. In particular we find using Eq. ([3]) 
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where T{z) is the Gamma function and where we used 

Tl(s)-1 

ps{o) - pm n ^D^^ (5) 

4=0 

n(s)-l 

Ps^Ds = Ps^Ds TT -Bd's- (6) 
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FIG. 2: (Color online) Power-law exponent zi vs. {A) accord- 
ing to analytic estimate Eq. pO|l (line) and from fitting the 
asymptotic decay of the numerically computed exact evolu- 
tion, Eq. ([2]), for a single realization of scaling factors (circles). 
The distribution for the scaling factor As (Bs) is chosen to be 
uniform on an interval centered around (A) {{B) = 0.3) with 
width 0.12 (0.1) and is independent from Bs (As). 



The products run over the sites D^s encountered on the 
direct path from s to the root s and n(s) is the number 
of levels ascended. 

As we are interested in the asymptotic behavior of Pi (t) 
for i ^ oo, we need to find the right boundary of the strip 
of convergence of Eq.(|3]). The Gamma function has no 
poles with positive real value and thus we concentrate on 
the factor containing the sum over the sites s. Fixing z 
real, the infinite sum of products of random variables 
diverges with probability one (Til if f3{z) > 1/2, where 



diz) ^ min , , 



A 



(7) 



Using the constraints (i) and (ii) of our random variables 
A and B, one finds that P{z) is monotonically increasing 
from /3(z = 0) < 1/2 to /3(z ^ cx)) = 1. This implies 
that there is a critical zi with 



/3(zi) = 1/2. 



(8) 



Thus the survival probability asymptotically decays as 



Piit) ^ C{t)t- 



(9) 



with C{t) growing at most logarithmically [l6|. We stress 
that zi depends on the probability distributions of the 
scaling factors A and B only and not on the specific val- 
ues appearing in a given tree. The power-law exponent 
zi is thus realization independent. 

For not too wide distributions of A and B, e.g. 
((-glf)^) < 2 (-glf)^, one can show that the minimum 
in Eq.([7]) appears at tr = 1 and Eq. ([8]) simplifies to 



A X 



(10) 



i=0 



A numerical example is shown in Fig. [T] for random vari- 
ables A and B independent and uniformly distributed on 
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an interval in logarithmic scale with In^ G [—4.21, —3.17] 
and InB G [-2.02,-0.98], giving zi = 1.88. It clearly 
demonstrates that, even though initially quite different 
power laws emerge, asymptotically all realizations have 
the same power-law exponent. In addition, one observes 
that the averaged curve shows this asymptotic exponent 
beginning at much smaller times. 

Moreover, in this case of not too wide distributions, 
one can show that fluctuations around the asymptotic 
algebraic decay vanish for each realization, as observed 
in Fig. [TJa. This is due to an increasing number of sites s 
contributing significantly to Eq. ([3]) as time increases. We 
have numerically demonstrated this by finding that the 
inverse participation ratio ^s[ps{t)/ Pit)]"^ of the con- 
tributing densities decreases like a power law as a func- 
tion of time (not shown). It stresses the fact that the con- 
tribution to P(t) at large times arises from an increasing 
number of paths [13] • 

In Fig. [2] we show numerical calculations of the 
power-law exponent fitted from the exact evolution 
Eq. ([2]) for a single realization of scaling factors, showing 
good agreement with the analytical prediction, Eq. (jTU]). 
This agreement gives also an indirect confirmation of 
the validity of the approximation used in Eq. Q. We 
note, that Eq. (jlOp reproduces the power-law exponent 
of deterministically chosen scaling factors A and B in 
Ref. 0. 

We now want to argue why this model is able to cap- 
ture the essential features of transport in Hamiltonian 
systems with a mixed phase space. The discrete nature of 
the Markov model is justified by the presence of invariant 
structures in phase space that act as partial barriers (e.g. 
cantori with minimum flux) and that effectively separate 
regions where the trajectory spends a long time before 
moving to the adjacent one. The specific set of transi- 
tion probabilities depends on the precise structure of the 
phase space, e.g., the transition probability for a trajec- 
tory to leave a region of area a leaking through a partial 
barrier with flux Aw is given by Aw /a. Area conserva- 
tion assures that the flux is the same in both directions 
and thus the scaling factors As in Eq. H]) are related to 
the scaling of areas of the connected regions, while the 
factors Bs are related to the ratio of escape rates. 

A typical mixed phase space is hierarchically orga- 
nized [171 ]. While approaching the boundary circle of a 
given island of regular motion (level-renormalization [l8| ) 
there are secondary island chains encircling the main one 
(class- renormalization [l^). This picture repeats on all 
scales ad inflnitum. Using these features of a mixed phase 
space Meiss and Ott modeled the transport in Hamil- 
tonian systems as a random walk on a binary tree with 
fixed scaling factors for level- and class-scaling, respec- 
tively. However, typical scaling factors in class-scaling 
and, at least initially, in level-scaling do vary and a more 
realistic model must take into account these variations. 



as it was suggested in Ref. 0- 

Our proposed model with random scaling factors, 
Eq. ([1]), is the simplest of this kind. It is not aimed 
at giving a quantitative description, but it elucidates the 
qualitative consequences of randomness in the scaling fac- 
tors, namely the presence of a well-deflned asymptotic al- 
gebraic decay with a realization-independent power-law 
exponent. We do not expect qualitative changes of this 
scenario by further refinements of the model, like partly 
fixing scaling factors for the modelin g of level-scaling or 
by changing the number of branches [l5| . 

In order to draw conclusions from our model about 
Hamiltonian systems, we make use of the hypothesis that 
there exists a universal mechanism generating the hierar- 
chical fine-scale structure of Hamiltonian systems, lead- 
ing to a universal distribution of scaling factors. In par- 
ticular, we consider different Hamiltonian systems, with 
respect to their fine-scale structure, as different realiza- 
tions from one statistical ensemble, just as this is the case 
for a specific realization of the random scaling factors A 
and B in the tree model. The immediate consequence of 
this hypothesis and of the realization independence of the 
power-law decays found in our model is, that all Hamil- 
tonian systems should exhibit the same universal power- 
law decay. The numerically observed different power-law 
decays over finite times are merely realization-dependent 
fiuctuations around the universal decay. 

As another important consequence, we are now able to 
effectively address the practical question of how to nu- 
merically determine the universal power-law exponent. 
The straightforward approach is to extract it from P{t) 
of a single Hamiltonian system determined up to very 
large times. However, in order to have a reliable power- 
law exponent from a single P{t) one would need much 
larger times than in previous studies (which would go 
beyond current computational resources). The reason 
is that the trajectories of the chosen Hamiltonian sys- 
tem must have explored a statistically significant sample 
of the fine structure of its phase-space in order to sam- 
ple the universal distribution of scaling factors. Instead, 
we propose to average P{t) over different Hamiltonian 
systems. This procedure averages out fluctuations and 
should give a much better defined power-law decay al- 
ready for smaller times, see e.g. Refs. [l|, [2^. It is im- 
portant to stress, that due to the results of our model 
such an averaging has a clear meaning: It is not an av- 
erage over different power laws, but it corresponds to a 
faster reduction of the realization-dependent fluctuations 
around the universal decay, obtained by an increased sta- 
tistical sampling of the universal distribution of scaling 
factors. 

Guided by this idea we study an ensemble of area pre- 
serving maps of the form pn+i = Pn - V'{qn); Qn+i = 
Qn +Pn+i with V{q) = §(7^ + for 100 pairs of param- 
eters (1 < a, 6 < 2). These maps have a stable fixed point 
at (0,0) and an unstable fixed point at (—^,0). The fine 
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FIG. 3: (Color online) (a) Survival probabilities P{t) for 100 
area preserving maps (thin lines), defined in the text, and 
their logarithmic average (thick line). Inset: phase space of 
one of these maps, showing the main island and the border 
line to the escape region, (b) Rescaled survival probabilities 
P{t)t'^ with the fitted exponent z = 1.57 from (a), (c) Vari- 
ance cr^(t) of \n[P{t)t'\ showing saturation. 



scale structure at the border of the regular island around 
the stable fixed point is quite different for all realiza- 
tions. We start 10^^ trajectories uniformly in a box next 
to the unstable fixed point, q G [— f , — | -I- 5],p S [— <5] 
with 5 = 0.01, and choose the escape region q < —a/b. 
Fig.Oshows the individual survival probabilities P{t) to- 
gether with their average. Note, that the invariant initial 
condition in the tree model corresponds in the Hamilto- 
nian system to a uniform initial distribution just in the 
chaotic component. Numerically, this initial condition is 
not achievable. Instead, we choose an initial condition 
in a box in the chaotic region far away from the island, 
which changes the exponent by one, z 



21|. 



zi + 1 

The averaged curve in Fig. [3] presents a well de- 
fined power-law decay with exponent z ~ 1.57 ± 0.03, 
where the error describes the variation for different fit- 
ting ranges. Quite importantly, the variance of the fluc- 
tuations around this averaged power-law saturates. This 
agrees with the prediction of our model (see Fig.[TJ:) and 
it gives direct evidence that all individual power-law de- 
cays are consistent with a single universal decay. We 
find similar results with the same exponent for an en- 
semble of kicked double-harmonic potentials. Motivated 
by our work, Altmann studied an ensemble of mod- 
ified standard maps, confirming our results qualitatively 
and quantitatively. 

In conclusion, we conjecture a universal power-law de- 
cay P{t) ^ from the analysis of a Markov model 
with random scaling factors and from the hypothesis of 
universality in the statistical properties of the fine-scale 
structure of typical Hamiltonian systems. Our numeri- 



cal investigations support this hypothesis and suggest the 
universal exponent z ~ 1.57 ± 0.03. Future studies will 
concentrate on a direct extraction of the universal distri- 
butions of area and escape-rate scalings in generic Hamil- 
tonian systems (both numerically and analytically) , that 
will eventually lead to a fully theoretical derivation of the 
universal exponent z. 

We thank E.G. Altmann for useful discussions. This re- 
search was supported in part by the DFG under contract 
KE537/3-3 and the NSF under Grant No. PHY05-51164. 
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